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It is shown that Guoping Zhang's results [G.P. Zhang, Phys. Rev. B 68, 153101 (2003)] for the 
charge-density-wave phase boundary in the half-filled one-dimensional extended Hubbard model are 
incorrect and that his criticism of my work [E. Jeckelmann, Phys. Rev. Lett. 89, 236401 (2002)] is 
groundless. 
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In Ref. Guoping Zhang presented density-matrix 
renormalization group (DMRG) results which contradict 
my DMRG calculations^ and Hirsch's quantum Monte 
Carlo (QMC) simulations^ for the charge-density-wave 
(CDW) phase boundary in the one-dimensional extended 
Hubbard model at half filling. In this Comment I show 
that Guoping Zhang's results are inaccurate and that his 
criticism of my work is groundless. 

Although the phase diagram of the extended Hubbard 
model is still partially controversial (see Refs. VJWmH and 
references therein), the CDW phase boundary V C (U) in 
the parameter space (U, V) was determined years ago&Si 
and has not been disputed in recent studies^*^* 7 ^ In 
Fig. ^ I show the results of various numerical investiga- 
tions for V C {U)—U /2 in the weak to intermediate coupling 
regime. There is an excellent overall agreement between 
Hirsch's QMC simulations,- the exact diagonalizations 
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FIG. 1: (Color online) Results for the CDW phase bound- 
ary V C (U): QMC simulations (Ref. |3j) (right triangle), exact 
diagonalizations (Ref. 6) (left triangle), level crossing anal- 
ysis (Ref.0) (circle), SSE-QMC simulations (Refs. H and H 
(up triangle), author's DMRG calculations [from Ref . (di- 
amond) and new results (square)], Yuzhong Zhang's DMRG 
calculations (Ref.0) (star), and Guoping Zhang's DMRG cal- 
culations (infinite-system algorithm) in Ref.Q (open down tri- 
angle) and in a previous work (Ref.|9|) (solid triangle). Dashed 
lines are guides for the eye. 
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FIG. 2: (Color online) Same results as in Fig. 1 but displayed 
using Zhang's representation (U/V c vs U). 



of Cannon et air Nakamura's level crossing analysis, 7 
the stochastic series expansion QMC (SSE-QMC) simu- 
lations of Sandvik et <d.J*- Yuzhong Zhang's DMRG cal- 
culations^, and my DMRG calculations 2 . In particular, 
my results agree quantitatively^ with the most recent 
and accurate numerical simulations^*^* 7 ^ Only Guoping 
Zhang's DMRG data^ deviate systematically from the 
other results. Therefore, there is clearly a problem with 
his calculations. 

The discrepancy between the various DMRG calcula- 
tion s] -^i 5 : 9 is not surprising. Guoping Zhang uses the 
infinite- system DMRG algorithm while Yuzhong Zhang 
and I use the more accurate finite-system DMRG al- 
gorithm. 11 It is well knowni^ that for many problems 
the infinite-system algorithm yield incorrect results while 
the finite-system algorithm gives essentially exact (nu- 
merical) results. In particular, it is essential to use the 
more reliable finite-system DMRG algorithm for inhomo- 
geneous systems such as a CDW ground state. Therefore, 
the discrepancy between Guoping Zhang's results and all 
other works just demonstrates the failure of the standard 
infinite-system DMRG algorithm for the present prob- 
lem. (See Ref. for another example of the infinite- 
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system algorithm failure and Ref. ^| for the successful 
investigation of the same problem with the finite-system 
algorithm.) 

In his paper Guoping Zhang wrongly claimed that 
my DMRG calculations (and the QMC simulations of 
Ref.0) failed to reproduce the weak coupling limit result 
V C (U) = U /2. In Refs. 00 the investigation of the phase 
diagram was focused on the intermediate- and strong- 
coupling regimes (i.e., U > 2t) and no analysis of the 
weak-coupling limit U <C t was performed. Here I present 
additional results for V C (U) calculated with DMRG for 
weaker couplings: V c /t = 0.260 ± 0.003 for U/t = 0.5, 
V c /t = 0.545±0.005for U/t = 1, and V c /t = 0.835±0.005 
for U — 1.5i. Moreover, I have calculated V C (U = 2t)/t 
more accurately and found 1.113 ± 0.005 (in agreement 
within the error bars with the value given in Ref. 0). 
These results are shown in Fig. ^ as square. They 
agree perfectly with other works^i and, clearly, they 
approach the weak-coupling result V C (U) = U/2 in the 



limit C/ — ^ 0. In Figs. 1 and 2 of Ref. Q] Guoping Zhang 
used a different representation of the data, U/V C (U) vs 
U, to analyze the weak-coupling limit. In Fig.0 I show 
again all data of Fig. fusing this representation. Clearly, 
the minimum of the ratio U/V C (U) occurs for U slightly 
smaller than 2t and the weak-coupling limit is recovered 
only for U smaller than 2t. Therefore, the U — ► limit of 
U/V C (U) cannot be determined using numerical data for 
U > 2t and the Figs. 1 and 2 of Ref. Q axe misleading. The 
DMRG and QMC data for U > 2t presented in Refs. @ 
and 3 arc fully compatible with the weak-coupling limit 
Vc{U) = U/2, contrary to Guoping Zhang's assertion in 
Ref. m 

In summary, comparisons with the results available 
in the literature confirm the accuracy and reliability of 
the DMRG calculations presented in Ref. |j. Guoping 
Zhang's results and conclusion are faulty due to the in- 
appropriate use of the infinite-system DMRG algorithm. 
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